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A Brownian Motor is a nanoscale or molecular device that combines the effects of thermal noise, 
spatial or temporal asymmetry, and directionless input energy to drive directed motion. Because of 
the input energy, Brownian motors function away from thermodynamic equilibrium and concepts 
such as linear response theory, fluctuation dissipation relations, and detailed balance do not apply. 
The generalized fluctuation-dissipation relation, however, states that even under strongly thermody- 
naniically non-equilibrium conditions the ratio of the probability of a transition to the probability 
of the time-reverse of that transition is the exponential of the change in the internal energy of the 
system due to the transition. Here, we derive an extension of the generalized fluctuation dissipation 
theorem for a Brownian motor for the ratio between the probability for the motor to take a forward 
step and the probability to take a backward step. 
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A Brownian Motor is a nanoscale or molecular device 
that combines the effects of thermal noise, spatial or tem- 
poral asymmetry, and directionless input energy to drive 
directed motion [J, 0, S] ■ Many biological motile systems 
may be driven by Brownian motors [j| , and chemists have 
been able to synthesize molecules that function as Brow- 
nian motors 0, 01 • In solution, viscous drag and thermal 
noise dominate the inertial forces that drive macroscopic 
machines. Because of the strong viscous drag, the mo- 
tion of such a Brownian motor is over-damped and in one 
dimension can be described by the simple equation 0] 



Ra-X = e{t) 



(1) 



where e{t) is Gaussian noise with mean fj, = and vari- 
ance cr^ = 2RkBT/dt, and R is the coefficient of viscous 
friction. In the following we use units where the thermal 
energy /cbT = 1- The generalized force X = X{a,'i(j{t)) 
can be written as the gradient of a scalar potential 
X = —dH/da where 



H{a,^{t)) = U{a)+ijj{t)z{a) 



(2) 



is the sum of an intrinsic potential due to chemical in- 
teractions and any external load and an external time 
dependent forcing term that is the product of canoni- 
cally conjugate intensive and extensive thermodynamic 
parameters z(a) and tpit), respectively The con- 

jugate parameters include, e.g., molecular volume and 
pressure, entropy and temperature, or dipole moment 
and field. The underlying system is typically spatially 
periodic (possibly with a homogeneous force or load F) 
so that U{a + L) = U{a) + AU, where AU = FL, and 
z{a + L) = z{a). 

For any fixed value of ip detailed balance requires 



P(ai + L,T\ 



PHa,,T\---\a,+L,Q) 



(3) 



where P{ai+L, T\ ■ ■ ■ \ai,0) is the conditional probability 
density that a particle starting at position ai at time 
goes to position ai+L at time T by the specific trajectory 
(sequence of positions and times) denoted by • • • , and 
P^ai,T\ ■ ■ - lai + L,0) is the conditional probability to 
follow the reverse of that process. The ratio depends only 
on the difference in energy between the initial and final 
points. It further holds that 



PiO,T\L,0) 



-AU 



(4) 



where the net probability P{L, T|0, 0) = Jq - ■ ■ P{oii + 
L,T| • • • |ai,0) is the integral over all trajectories from 

(0,0) to (L,r). 

A time dependent modulation, ijiit), causes dissipation 
and breaks detailed balance, in which case Eqs. (3) and 
(4) do not hold. It is even possible to have 



P{L,T\Q,0) 
P(0,r|L,0) 



> 1; 



e-'^^ < 1 



(5) 



where the external stimulus ^{t) provides energy to drive 
uphill motion [^, [l3| . 

The generalized fluctuation-dissipation theorem [ill . 
12| states that even under strongly thermodynamically 
non-equilibrium conditions the ratio of the probability of 
a forward (F) transition to the probability of the time- 
reverse (R) of that transition is the exponential of the 
change in the internal energy of the system due to the 
transition 



Pf(J^,T|...|0,0) 

PR(o,r|...|L,o) 



W-AU 



(6) 



where W is the work supplied to the system by the ex- 
ternal modulation in the forward trajectory. Here, an 
extension of the generalized fluctuation dissipation theo- 
rem for a Brownian motor is derived to obtain the ratio 
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FIG. 1: Depiction of symmetry related trajectories of a Brow- 
nian particle in a periodic ratchet potential, a) Snapshot 
of the potential described in Eq. (22) with a particle at 
a = 0. For any constant ip, the forward (F) and backward 
reverse (BR) trajectories are identical (and Pp = Pbr), as 
are the time reverse (R) and backward (B) trajectories (and 
Pr = Pr) and by detailed balance Pp/Pb ~ e~ AU and the 
net motion of the particle is to the left. External modula- 
tion tp{t) breaks detailed balance and can drive net motion 
to the right against the load, b) With an external forcing 
^{t) all four trajectories - F,R,B, BR - are distinct from one 
another. As described in the text, an extension of the general- 
ized fluctuation dissipation relation can be used to derive the 
ratios of the probability densities for these symmetry related 
trajectories. 



generalized force, Xj+i = Xf_|_j + X°_^^, are 

where Af^_^_l = /'^(ai+i) - f^{ai) for / = U, z and 
k = e,o. For every forward trajectory {a{t),'ip{t)} with 
probabihty Pp, defined by 



i< = > ai > 



L, 



M-l 



-Pf = P(ai+i|ai;?/'i+i), 



(9) 



i=0 



there are three symmetry related trajectories. One is a 
time reverse trajectory {a{—t),^{--t)} obtained by 
switching the sign of time. For a time periodic system 
reversing time is equivalent to the transformation t 

i-r-t) , 



ti = L > Qfin-l > 



M-l 



— > ai — > 0, 



pR = n P{ai\ai+i;ipi+i) 



(10) 



i=0 



Another is a backward trajectory {—a{t), ip{t)} obtained 
by switching the sign of the position variable. For a space 
periodic system this is equivalent to the transformation 

a ^ (L — a), 



B = 



between the probabihty for the motor to take a forward 
step and the probability to take a backward step in for- 
ward time. 

First, we write equation (1) as a more rigorous finite 
difference or update equation and convert to unit normal 
Gaussian noise A^(0, 1) [isj 



tti+i - a; - R-^Xi+iAt = V2R-^At N{0, 1) (7) 

The time interval At is chosen to be sufficiently short 
that the change in position Aa is very small. We used 
the relation N{^i,a^) = n + aN{0,l) where A^(0, 1) is 
a Gaussian random variable with zero mean and unit 
variance the values, n, of which occur with probability, 
P{n) = exp (— n^/2)/(-\/27r). Any two values of n are 
uncorrclatcd (nirik) = (5i^k- 

Broken symmetry is an essential feature of a Brownian 
motor, so we split each of the position dependent terms 
into even and odd components U{a) = U^{ct) -I- U°{oi) 
and z(a) = z'^ipi) + 2° (a), where for any function 
p{-a) = f(a) and /°(-a) = -/°(a). Finite differ- 
ence expressions for the even and odd components of the 



M-l 

-Pb = n ^(-"i+il (11) 

i=0 

The third is a backward reverse trajectory 
{—«(<),■(/'(—*)} obtained by switching the sign of 
both time and of the position variable. For a time 
and space periodic system this is equivalent to the 
transformation a — > (i — a) and t — > (T — t), 

T3T3 — T '''m V'm-l ■4'2 lAl r, 

BK = — P — > — Q'in-i — * ■ ■ ■ — > — ai — * 0, 
M-l 

Pbr = n ^(""il " "i+i; ^i+i): (12) 

i=0 

Viewing Eq. (7) as a rnap ping between the "noise" 
space and "position" space [l5|, ll6| , the conditional prob- 
ability density given that the system is at position ai after 
the i*'' step, and that the value of the field is for the 
(i+l)-'** step is seen to be 

-(Aa-J^^^X^^j^At)'^ 

P(ai+i|ai,V'i+i) = . „ • (13) 
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where Aa = (ofi+i — ai). The ratio of the probability 
density for the forward and time reverse step is 



P(Q;i+i|o;i,V>i+i) 
P(ai|ai+i,'i/'i+i) 



(14) 



and the ratio between the forward and time reverse tra- 
jectory is 



Pf 



where 



Wf 



M-l 

E 

i=0 



(Az^(a) + Az°(a)) 



(15) 



(16) 



is the total external work done in the forward trajec- 
tory. Eg (15) is the generalized fluctuation dissipation 
relation ll|, [iJI and the change in the internal energy 
of the system, AE = AU — Wp, is the dissipated work. 
The ratio of the probability density for a backward and 
backward time reverse is similarly obtained, 



Pb 
Pbk 



where 



M-l 



Wb=Y^ Vi+i {Az%a) - Az°(a)) 



(17) 



(18) 



i=0 



is the total external work done in the backward trajec- 
tory. Finally, the ratio between a forward and backward 
step is 



P(Qi+l|ai,^i+l) _ X'+iAa-X.<ViA-°+ii?-iAi 



P(-Q!i+l| - ai,l/'i+i) 



(19) 



The central result of this paper, the ratio for the prob- 
ability densities for a forward and backward trajectory, 
follows immediately 



P^ 
Pb 



-Br^ X°X°dt 



(20) 



where we have taken the limit At — > to get the inte- 
gral form. Equation (20) highlights the importance of 
broken symmetry - if either X° or X° is zero, the ratio 
of the probability for a forward step to a backward step 
is governed solely by the homogeneous force acting on 
the system and is independent of the work pumped in by 
the time dependent modulation. Consider the standard 
ratchet potential below Q 



U{a) = Uo cos (2a) -f Fa 



z{a) = zq [cos (stt) cos (a) + sin (stt) sin [a)] 



(21) 



where the asymmetry parameter is — 1 < s < 1. For 
Uq > tpo^o^ where "00 is the amplitude of the external 
modulation, a single spatial period of the potential has 
two clearly defined energy wells, say A and B, where the 
sign oiijj determines the relative energies of the two wells, 
and the relative heights of the barriers between them. We 
can then describe the motion as a random walk between 
the two wells, ^ ^ , where a clockwise transition 

indicates a half-step to the right, and a counterclockwise 
transition indicates a half-step to the left. The transition 
constants are 

fcAB = fco [e-^^(<^'=0°)]'/' 
lBA = fco[e-^^/(0^(/>°)]'^' 
fcAB^fco [e^^(0V0°)]'^^ 



(22) 



where cj)" = e4^»'>(*) (''^) and 0° = e'*^'"^(*) '^''^ . Irre- 
spective of the value of s or the form of ip{t) a corollary 
of detailed balance for rate processes. 



'^AB'tBA _ g-AC/ 



SABKBA 



(23) 



holds at every instant. The net motion can be solved 
analytically for small amplitude [l7| '4'it)^ ^-nd in spe- 
cial cases, such as square wave perturbation, for arbi- 
trary amplitude 3 '^(i)- ^ similar model has been 
proposed for adiabatic transport [l3|- For the specific 
case that 'ip{t) is externally generated dichotomic noise 



_v]/ ) in which i/j{t) switches between +4' and 



— 5* with a Poisson distributed random lifetime (average 
1 /"/) ( a situation particularly relevant for Brownian mo- 
tors that are driven, e.g., by the stochastic binding of 
chemical fuel molecule and release of product) the com- 
bined stepping/switching process can be described by a 
single diagram Q 



B,+ 



B,- 



(24) 



This picture emphasizes the idea of a minimal Brow- 
nian motor as two coupled two-state processes. One 
process is the externally driven dichotomic modulation 

7 

j^q, ^ and the other is the thermally activated 

7 

Stepping ^ ^ . 

The overall diagram can be broken into six cycles Q 
- two cycles for the uncoupled stepping, one with fixed 
+5* and the other with fixed —4', two cycles for the dis- 
sipative back and forth motion with no net stepping to 



4 



the right or left, and two cycles describing net stepping 
coupled to the external fluctuation. The last two, cou- 
pled, cycles are of most interest. The forward, reverse, 
backward, and backward reverse paths are 



B, 



B, 



B,- 



A,+ 

A,+ 
A,- 



B, 



B,- 



(25) 



BR 



B. 



B,- 



The probability for completion of a cycle is proportional 
to the product of the transition constants in the cycle Q . 
The proportionality constants involve rate constants for 
back and forth transitions, lifetimes of the states within 
the cycle, etc. Importantly, since F, R, B, and BR direc- 
tional cycles (Eq. (26)) involve the same states, the pro- 
portionality constants arc the same for all of these sym- 
metry related cycles. Thus, with h^j^ = kij{ip{t) = i^P), 
it is easy to derive 



Pp ^ kA+B+kB-A- _ ^W-AU 



(26) 



Pb _ kA+B+kB-A- 



^W+AU 



^BR ks.A^kA^B. 



where IF = In [4)°^.) is the work done in the forward cycle 
when the energy is increased by 2zo^' sin (stt) in going 

from A, > ^, + and again from B, H > B,— . Here 

and below = (/)('='°)(i/)(t) = +*). The ratio of the 

probabilities for a forward and backward cycle is 



Pf 
Pb 



kA+B+ ks^A- 
kA+B+ ks^A- 



-AUic 



(27) 



and the ratio of the net forward to backward steps is 



Pf + Pbr 
Pb + Pr 



-AU 



1 



(28) 



The expansion of the coefficient in Eq. (29) involves only 
even powers of the amplitude ipo of the external driv- 
ing - the Brownian motor mechanism is a fundamentally 
non-linear effect of the external driving [3l ■ It does 
not necessarily require a large amplitude driving to ob- 
serve experimentally (20| however since the linear term 
is, by symmetry, identically zero. The non-monotonic 
frequency response observed experimentally [20| and ex- 
plained theoretically 2l| arises from the uncoupled tra- 



jectories which take on greater or lesser importance de- 
pending on the frequency of the applied signal. The sym- 
metry relations derived here are evidently frequency in- 
dependent. 



Many recent synthetic implementations of molecular 
brownian motors involve motion between discrete bind- 
ing sites. Without an external driving the thermally acti- 
vated transitions show no long time order irrespective of 
structural asymmetry, in consistency with the principle 
of detailed balance. By using external energy to ma- 
nipulate the environment, even in a seemingly random 
way, it is possible to break detailed balance and to drive 
directed motion. The symmetry relations derived here 
remain valid even under the action of an external pertur- 
bation and provide insight into how it may be possible 
to optimize synthetic molecular motors. 
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